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LONG-TERM  GOALS 

To  determine  the  mRNA  and  protein  expression  patterns  of  host  macrophages  during  anthrax  infection. 

OBJECTIVES 

Objectives  are  to  develop  the  technologies  to  do  high-level  functional  genomics,  proteomics  and 
bioinformatics  of  host  macrophages  during  early,  establishment  stages  of  anthrax.  The  approaches  are 
designed  to  track  each  macrophage  gene's  activities  within  the  entire  host  genome,  without  prejudice. 
In  this  manner  it  is  anticipated  that  differentially  expressed  host  genes  may  play  a  direct  role  in  the 
pathogenesis  and  virulence  of  the  disease  anthrax.  The  complete  determination  of  proteins  constituting 
the  host  proteome  during  anthrax,  with  emphasis  on  immune  cells,  will  be  performed  and  used  to 
determine  anthrax-specific  immune  responses  and,  hopefully,  key  break  points  in  anthrax-related 
pathologies. 

APPROACH 

To  perform  host  cell  functional  genomics  and  proteomics  we  exploit  the  human  and  mouse  genome 
sequences.  Total  RNA  and  protein  is  collected  from  host  cells  (control)  and  B.  anthracis- challenged 
host  cells  (experimental)  as  the  relevant  culture  model  for  early  anthrax.  Protein  and  mRNA  species 
are  identified  and,  when  possible,  quantified,  by  microarray  and  LC-LC-tandem-mass  spectroscopy 
(MudPIT)  methodologies.  This  is  a  multidisciplinary  effort  with  the  experimental  “wet”  studies  being 
performed  in  my  laboratory,  the  microarray  (gene-chip)  RNA  analyses  being  performed  at  TIGR  and 
the  University  of  Michigan  and  the  tandem  mass  spectroscopy  protein  determination  being  performed 
atjScripps;  with  all  3  groups  performing  bioinformatics  analyses.  The  large  amounts  of  data  generated 
represent  the  entirety  of  all  host  genes  and  proteins  as  determined  from  the  genome  sequences.  Data  is 
analysed  by  available  sifting  and  grouping  software  programs  and  those  RNA  and  protein  species 
showing  differential  expression  patterns  related  to  pathogenesis  are  studies,  individually,  using  current 
molecular  pathogenesis  technologies. 

WORK  COMPLETED 

Complete,  validated  host  mRNA  and  total  protein  expression  of  macrophage  data  sets  from  0  hour,  3 
hour  and  6  hour  post  anthrax  infections  and  post-anthrax  lethal  toxin  challenge  (and  controls)  were 
collected  with  high  confidence,  compared  and  analysed.  Gene  expression  data  sets  were  published 
(attached)  in  2005  [Bergman,  et  al.  Murine  Macrophage  Transcriptional  Responses  to  Bacillus 
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anthracis  Infection  and  Intoxication,  Infection  &  Immunity,  73:1069-1079.].  Complete,  curated 
mRNA  data  sets  were  deposited  at  the  TIGR  website  (htpp://pga.tigr.org),  the  ArrayExpress  data  base 
(accession  number  E-TIGR-84)  and  submitted  as  a  public  resource  to  the  NIAID  Biodefense 
Proteomics  Research  Centers.  Proteomics  data  sets  are  submitted  to  the  NIAID  Biodefense  PRCs. 

RESULTS 

Findings  determined  3  categories  of  macrophage  mRNA  expression  patterns. 

1 .  Macrophage  mRNA  expression  unchanged  during  anthrax  infections 

2.  Macrophage  mRNA  expression  changes  during  anthrax  infections  that  are  similar  to  what  is 
known  by  infection  with  most  other  pathogens. 

3.  Macrophage  mRNA  expression  changes  during  anthrax  infections  that  are  unique  to  anthrax. 

4.  These  data  sets  were  made  available  to  the  scientific  community  as  a  general  resource. 

Several  hundreds  of  host  genes  were  identified  that  were  differentially  expressed  during  the  early 
stages  in  the  model  system  for  anthrax.  These  loci  included  both  genes  that  are  known  to  be  regulated 
differentially  in  response  to  many  other  bacterial  pathogens  as  well  as  those  that  appear  to  be 
differentially  regulated  in  response  to  B.  anthracis  but  not  other  bacterial  species  that  have  been  tested 
to  date.  Our  data  indicates  a  transcriptional  basis  for  a  variety  of  physiological  changes  observed 
during  infection,  including  the  induction  of  apoptosis  caused  by  the  infecting  bacteria.  The  expression 
patterns  underlying  B.  anthracis- induced  apoptosis  led  us  to  test  further  the  importance  of  one  very 
highly  induced  host  gene,  that  for  ornithine  decarboxylase  (ODCT  Results  indicate  that  this  enzyme 
plays  an  important  and  previously  unrecognised  role  in  suppressing  apoptosis  in  B.  anthracis- infected 
cells.  We  postulate  that  ODC  can  be  a  useful  diagnostic  marker  for  pre-symptomatic  anthrax 
infections. 

These  data  were  published: 

Bergman,  N.,  Passalacqua,  K.,  Gaspard,  R.,  Shetron-Rama,  L.,  Quackenbush,  J.,  &  P.  Hanna. 
2005.  Murine  Macrophage  Transcriptional  Responses  to  Bacillus  anthracis  Infection  and 
Intoxication.  Infection  &  Immunity.  73:1069-1079. 

Parallel  to  the  mRNA  data,  MudPIT  proteomics  data  have  collected  on  anthrax-challenged 
macrophages.  These  data  have  been  analysed  and  the  data  sets  submitted  to  the  NIAID’s  Biodefense 
Proteomics  Research  Centers  as  a  resource  for  the  scientific  community.  Follow-up,  validation, 
studies  continue  with  NIAID  support. 

IMPACT/APPLICATIONS 

The  host  responses  to  any  infectious  agent  can  dictate  microbial  clearance  or  disease  or,  in  the  case  of 
anthrax,  death.  The  use  of  anthrax  as  the  model  system  in  developing  these  technologies  has  many 
advantages.  It  is  considered  an  important  candidate  in  the  instance  of  BW  and  the  data  generated 
should  allow  real  improvements  in  rapid  detection  and  diagnostics  technologies,  especially  in  the  arena 
of  anthrax  spore-specific  detection  and  now,  potentially,  in  discovery  of  candidates  for  new 
diagnostics.  Discovery  of  several  key  intervention  points,  at  the  level  of  points  for  development  of 
immunomodulators,  are  also  indicated  by  our  data  but  require  yet  validation  (ongoing). 
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TRANSITIONS 

None  yet. 

RELATED  PROJECTS 

The  ONR  sponsored  3  separate  grants  on  this  topic  in  a  close  collaboration  between  our  lab  (anthrax 
pathogenesis),  TIGR  (anthrax  genome  sequencing  and  microarray  development)  and  Scripps 
(functional  proteomics).  Cross  training  and  an  interdisciplinary  approach  was  emphasised.  The 
projects  continue  with  current  support  from  the  NIAID  (through  ’09).  We  are  one  of  seven  NIAID 
Biodefense  Proteomics  Research  Centers. 
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Interactions  between  Bacillus  anthracis  and  host  macrophages  represent  critical  early  events  in  anthrax 
pathogenesis,  but  their  details  are  not  clearly  understood.  Here  we  report  the  first  genomewide  characteriza¬ 
tion  of  the  transcriptional  changes  within  macrophages  infected  with  B.  anthracis  and  the  identification  of 
several  hundred  host  genes  that  were  differentially  expressed  during  this  intracellular  stage  of  infection.  These 
loci  included  both  genes  that  are  known  to  be  regulated  differentially  in  response  to  many  other  bacterial 
pathogens  and  those  that  appear  to  be  differentially  regulated  in  response  to  B.  anthracis  but  not  other  bacterial 
species  that  have  been  tested.  These  data  provide  a  transcriptional  basis  for  a  variety  of  physiological  changes 
observed  during  infection,  including  the  induction  of  apoptosis  caused  by  the  infecting  bacteria.  The  expression 
patterns  underlying  B.  anthracis- induced  apoptosis  led  us  to  test  further  the  importance  of  one  very  highly 
induced  macrophage  gene,  that  for  ornithine  decarboxylase.  Our  data  show  that  this  enzyme  plays  an  impor¬ 
tant  and  previously  unrecognized  role  in  suppressing  apoptosis  in  B.  anthracis- infected  cells.  VVe  have  also 
characterized  the  transcriptional  response  to  anthrax  lethal  toxin  in  activated  macrophages  and  found  that, 
following  toxin  treatment,  many  of  the  host  inflammatory  response  pathways  are  dampened.  These  data 
provide  insights  into  B.  anthracis  pathogenesis  as  well  as  potential  leads  for  the  development  of  new  diagnostic 
and  therapeutic  options. 


Bacillus  anthracis  is  a  gram-positive,  endospore-forming  bac¬ 
terium  that  is  the  causative  agent  of  anthrax.  This  disease  is 
primarily  zoonotic,  affecting  grazing  herbivores,  and  naturally 
occurring  human  infections  are  rare  in  developed  countries 
(13).  Even  so,  the  resilience  of  the  B.  anthracis  endospore,  the 
potentially  lethal  nature  of  the  disease,  and  the  efficiency  with 
which  it  infects  via  an  aerosol  route  make  it  an  ideal  biological 
weapon,  and  for  these  reasons  the  organism  and  the  disease  it 
causes  have  recently  come  under  increased  scrutiny  (28,  29). 

Infection  by  B.  anthracis  begins  with  entry  of  dormant  en- 
dospores  into  the  body  through  one  of  three  possible  routes: 
cutaneous  exposure,  ingestion,  or  inhalation.  In  systemic  cases, 
spores  are  efficiently  scavenged  by  resident  phagocytes  and 
transported  to  the  regional  lymph  nodes  (16,  17,  24,  47).  The 
spores  germinate  en  route,  transforming  into  replicative  vege¬ 
tative  cells.  The  nascent  cells  multiply  within  the  phagocyte  and 
begin  to  produce  virulence  factors  such  as  the  two  binary 
exotoxins  lethal  toxin  and  edema  toxin,  eventually  killing  the 
host  cell  and  escaping  into  the  local  environment  (the  blood¬ 
stream,  in  systemic  anthrax).  Up  to  this  point,  the  symptoms  of 
anthrax  are  relatively  mild,  but  once  in  the  blood  the  bacteria 
grow  to  a  very  high  titer  (sometimes  approaching  108  organ¬ 
isms  per  ml  of  blood)  (13,  47),  and  this  combination  of  septi¬ 
cemia  and  toxemia  quickly  leads  to  shock  and  death  of  the 
host.  Because  anthrax  is  usually  diagnosed  at  this  late  stage, 
the  disease  is  very  difficult  to  control.  Understanding  the 
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events  prior  to  this  point — those  occurring  in  association  with 
the  host  phagocyte— is  therefore  critical  for  the  development 
of  new  therapeutic  options. 

Although  much  of  the  B.  anthracis -host  interaction  takes 
place  within  the  phagocyte,  the  individual  molecular  events 
that  make  up  this  interaction  are  not  well  understood.  Some  of 
the  factors  that  are  critical  for  growth  of  B.  anthracis  inside  the 
cell  are  beginning  to  be  understood  (23,  24,  30,  31),  but  the 
ways  in  which  the  macrophage  responds  to  the  growing  infec¬ 
tion  remain  unclear.  In  order  to  better  define  these  responses 
and  perhaps  identify  new  points  at  which  new  therapeutic  or 
diagnostic  options  might  be  possible,  we  sought  to  examine  this 
response  in  detail  on  a  global  scale  with  DNA  microarrays. 
This  approach  has  been  successful  in  several  recent  studies 
examining  the  response  of  a  variety  of  mammalian  cell  types  to 
a  range  of  infecting  bacteria,  including  both  intra-  and  extra¬ 
cellular  pathogens  (3,  4,  6,  38, 46).  These  analyses  have  yielded 
insights  into  the  common  immune  response,  and  they  have  also 
identified  genes  that  are  specifically  involved  in  the  response  to 
a  particular  bacterium. 

In  this  study  DNA  microarrays  were  used  to  examine  the 
transcriptional  response  of  the  murine  macrophage-like  RAW 
264.7  cell  line  during  infection  by  B.  anthracis  and  also  in 
activated  macrophages  following  treatment  with  B.  anthracis 
lethal  toxin.  In  comparing  the  expression  patterns  induced  in 
each  of  these  conditions,  we  anticipated  being  able  to  identify 
host  genes  specifically  induced  or  repressed  in  response  to 
toxin  exposure  as  well  as  genes  specifically  regulated  during 
infection.  Within  these  response  programs,  we  identified  genes 
that  are  part  of  the  common  immune  response  induced  by 
many  or  all  bacterial  pathogens  as  well  as  those  that  have  not 
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been  reported  previously  in  other  systems  and  thus  may  be 
unique  to  the  B.  anthracis- macrophage  interaction.  Our  find¬ 
ings  confirm  many  of  the  biochemical  and  physiological  phe¬ 
nomena  previously  observed  during  this  interaction,  including 
the  importance  of  apoptosis  in  the  pathogenic  process  and  the 
significant  modulation  of  cytokine  synthesis  by  B.  anthracis.  In 
addition,  we  have  identified  a  limited  number  of  macrophage 
response  genes  that  appear  to  be  different  from  the  common 
responses  to  previously  characterized  pathogenic  bacteria  and 
may  serve  as  useful  leads  in  the  development  of  future  diag¬ 
nostic  and/or  therapeutic  options  targeting  B.  anthracis  and  its 
close  relatives. 

MATER  LA  LS  AND  METHODS 

Growl h  of  bacteria  and  cell  culture  strains.  Bacillus  anthracis  Sterne  34 F2 
(containing  plasmid  pXOl  but  not  pX02)  was  grown  in  brain-heart  infusion 
broth  or  on  sheep's  blood  agar  plates,  and  spores  were  prepared  as  in  Liu  et  al. 
(36).  The  murine  macrophage-like  RAW  264.7  cell  line  (ATCC  TIB-71)  was 
maintained  in  Dulbecco’s  modified  Eagle’s  medium  (DMEM)  with  10%  fetal 
bovine  serum  (Gibco-BRL)  at  37°C  with  5%  C02  in  a  humidified  incubator. 
Cells  were  visualized  by  Diff-Quick  staining  (Baxter  Scientific). 

Infection  and  toxin  exposure  conditions.  Twenty-four  h  prior  to  infection, 
RAW  264.7  cells  were  counted  and  switched  from  DMEM- 10%  fetal  bovine 
serum  to  minimal  essential  medium  (MEM)  with  10%  horse  serum,  as  this 
combination  has  been  shown  to  inhibit  outgrowth  of  extracellular  spores.  Cells 
were  seeded  at  3  x  107  cells  per  30  ml  of  medium  in  150-cm2  cell  culture  flasks 
and  cultured  overnight.  Infections  were  begun  by  removing  the  medium  and 
adding  15  ml  of  fresh  MEM-10%  horse  serum  containing  B.  anthracis  Sterne 
34F2  spores  (a  strain  that  carries  pXOl  but  lacks  pX02  and  is  thus  able  to 
produce  both  binary  exotoxins  but  not  capsule)  at  a  multiplicity  of  infection  of 
10:1. 

In  experiments  in  which  toxin  exposure  was  being  examined,  the  procedure 
was  essentially  the  same  except  that  the  fresh  medium  contained  B.  anthracis 
protective  antigen  (500  ng/ml)  and  lethal  factor  (various  concentrations  for 
preliminary  experiments,  8  ng/ml  for  treatments  prior  to  RNA  isolation).  At 
these  concentrations  of  protective  antigen  and  lethal  factor,  previous  studies 
have  shown  that  essentially  all  of  the  lethal  factor  is  bound  by  protective  antigen, 
and  thus  the  functional  concentration  of  lethal  toxin  is  the  same  as  that  of  lethal 
factor  (5).  Note  also  that  free  protective  antigen  at  these  concentrations  did  not 
have  an  effect  on  cell  viability,  as  seen  in  our  own  control  experiments  (data  not 
shown)  and  previous  studies  (45).  Salmonella  lipopolysaccharide  (0.5  jAg/ml)  was 
added  simultaneously  with  lethal  toxin  when  noted.  For  control  experiments, 
fresh  medium  was  added  without  spores  or  toxin. 

RNA  isolation.  Macrophage  RNA  was  isolated  al  the  desired  time  points  by 
decanting  the  medium  and  extracting  sequentially  with  Trizol  (Invitrogen)  and 
bromochloropropane  according  to  the  manufacturer’s  instructions.  The  RNA 
was  then  precipitated  with  isopropanol  and  washed  with  75%  ethanol  and  then 
dried.  The  pellets  were  resuspended  in  200  p,l  of  H20.  The  RNA  was  further 
purified  with  RNeasy  Midi  columns  (Qiagen)  and  again  resuspended  in  H20. 
RNA  quality  and  purity  were  examined  by  measuring  absorbance  at  260  and  280 
nm  and  by  visualization  on  a  denaturing  agarose  gel. 

Generation  of  probes  for  microarray  experiments.  DNA  probes  for  microarray 
experiments  were  generated  as  described  by  Abulencia  et  al.  (I).  Briefly,  15  jxg 
of  total  RNA  was  reversed  transcribed  in  a  mixture  containing  6  p.g  of  random 
hexamers  (Invitrogen),  0.01  M  dilhiothreitol,  an  aminoallyl-deoxynucleoside 
triphosphate  mixture  containing  25  mM  each  dLATP,  dCTP,  and  dGTP,  15  mM 
dTTP,  and  10  mM  aminoallyl-dUTP  (Sigma),  reaction  buffer,  and  400  units  of 
Superscript  II  reverse  transcriptase  (Invitrogen)  at  42°C  overnight.  The  RNA 
template  was  then  hydrolyzed  by  adding  NaOH  and  EDTA  to  a  final  concen¬ 
tration  of  0.2  and  0.1  M,  respectively,  and  incubating  at  70°C  for  15  min.  Unin¬ 
corporated  aminoallyl-dUTP  was  removed  with  a  QIAquick  column  (Qiagen). 
The  probe  was  eluted  with  a  phosphate  elution  buffer  (4  mM  KP04,  pH  8.5,  in 
ultrapure  water),  dried,  and  resuspended  in  0.1  M  sodium  carbonate  buffer  (pH 
9.0).  To  couple  the  aminoallyl  with  fluorescent  labels,  normal  human  serum-Cy3 
or  normal  human  serum-Cy5  (Amersham)  was  added  in  the  dark  at  room 
temperature  for  I  h.  Uncoupled  label  was  removed  with  the  QIAquick  column 
procedure. 

Microarray  hybridization,  scanning,  and  image  analysis.  Ami  nosilane -coated 
slides  printed  with  a  set  of  27,648  target  sequences  were  prehybridized  in  5X  SSC 


(1  x  SSC  is  0.15  M  NaCI  plus  0.015  M  sodium  citrate)  (Invitrogen),  0.1%  sodium 
dodecyt  sulfate,  and  1%  bovine  serum  albumin  at  42°C  for  60  min.  The  slides 
were  then  washed  at  room  temperature  with  distilled  water,  dipped  in  isopro¬ 
panol,  and  allowed  to  dry.  Equal  volumes  of  the  appropriate  Cy3-  and  Cy5- 
labeled  probes  were  combined  and  dried  and  then  resuspended  in  a  solution  of 
50%  formamide,  5x  SSC,  and  0.1%  sodium  dodecyl  sulfate.  The  RNA  was 
supplemented  with  20  gg  of  mouse  COTl  DNA  and  20  p.g  of  poly(A)  DNA  and 
then  heated  (95°C,  3  min)  and  quickly  cooled  on  ice  for  30  s.  The  probe  mixture 
was  then  added  to  the  microarray  slide  and  allowed  to  hybridize  overnight  at 
42°C. 

Hybridized  slides  were  washed  sequentially  in  solutions  of  lx  SSC-0.2%  SDS, 
O.lx  SSC-0.2%  SDS,  and  0. lx  SSC  at  room  temperature,  then  dried  in  air,  and 
scanned  with  an  Axon  GenePix  4000  scanner.  Individual  TIFF  images  from  each 
channel  were  analyzed  with  TIGR  Spotfinder  (available  at  http://www. tigr.org 
/software/tm4)  (50).  Spots  whose  shape  deviated  significantly  from  a  true  circle 
and  those  with  signal-to-noise  ratios  of  <3.0  were  eliminated  from  further  anal¬ 
ysis. 

Microarray  normalization  and  data  analysis.  Microarray  data  sets  were  nor¬ 
malized  with  a  locally  weighted  linear  regression  (LOWESS)  algorithm  as  im¬ 
plemented  in  the  TIGR  MIDAS  software  package  (50).  This  program  also 
allowed  us  to  filter  the  data  in  order  to  remove  spots  for  which  flip  dye  replicates 
deviated  considerably  from  the  expected  reciprocal  relationship.  The  resulting 
data  were  then  log2  transformed,  and  adjusted  relative  to  the  time  zero  point 
(uninfected)  as  described  in  Baldwin  et  al.  (3).  In  brief,  this  was  done  by  sub¬ 
tracting  the  log  (2)  of  the  time  zero  point  from  each  of  the  other  data  sets,  such 
that  the  reference  RNA  was  effectively  eliminated  and  changes  in  gene  expres¬ 
sion  during  the  course  of  the  6-h  experiment  could  be  measured  relative  to  the 
zero  time  point.  Following  this  step,  the  data  were  averaged  over  duplicate  genes 
within  each  array  (when  present)  and  over  duplicate  arrays  done  for  each  ex¬ 
periment. 

Differentially  expressed  genes  were  determined  at  the  99%  confidence  level 
for  each  experimental  data  set  by  assuming  that  the  log2  ratios  for  each  data  set 
approximate  a  normal  distribution  (initial  characterizations  of  the  data  sets 
obtained  in  this  study  showed  that  the  data  did  not  show  an  appreciably  intensity- 
dependent  structure  after  LOWESS  normalization;  see  supplemental  Fig.  SI  and 
all  other  supplemental  figures  and  tables  mentioned  in  this  article  at  http: 
//www.sitemaker. umich.edu/hanna_Iab),  and  selecting  genes  with  log2  ratio  val¬ 
ues  that  deviated  by  >2.58  standard  deviations  from  the  mean.  In  almost  evei^ 
case  this  level  of  confidence  translated  to  a  >2-foId  change  in  expression  levels, 
but  in  the  rare  instance  that  this  was  not  true,  our  filtration  was  adjusted  such 
that  significant  differential  expression  was  defined  as  deviating  by  >2.58  standard 
deviations  and  changing  >2-fold  in  expression  level.  Filtration  for  significant 
differential  expression  and  further  analyses  were  done  with  MS  Excel. 

Gene  ontology  (GO)  annotation  of  genes  identified  as  differentially  expressed 
was  carried  out  with  the  Onlo-Express  software  package  (available  at  hup: 
//vortex.cs.wayne.edu/Projects.html)  and  the  UniGene  cluster  identifiers  for  each 
target  sequence  on  the  microarray.  Probabilities  of  over  represent  at  ion  were 
estimated  with  a  chi-square-based  probability  model,  with  Fisher’s  exact  test 
substituted  when  necessary  because  of  small  sample  size.  When  both  a  category 
and  its  subcategories  showed  statistically  significant  overrepresentalion,  the 
smallest  (most  specific)  functional  group  that  showed  significant  overrepresen¬ 
tation  was  reported, 

RT-PCR  assays.  Reverse  transcription  (RT)  and  PCR  were  used  to  validate 
the  microarray  data  corresponding  to  the  expression  profiles  of  selected  genes. 
Primers  were  synthesized  by  Invitrogen,  Inc.,  and  Integrated  DNA  Technologies. 
The  sequences  used  were  5'-ACATGATGCTGTAAACTTTATTTGCGA  and 
5'-CCTCTTTTGTTCCGCTATCGATC  (methylmalonate-semialdehyde  dehy¬ 
drogenase);  5  -CAGCAGGCTTCTCTTGGAAC  and  5'-CATGCATTTCAGG 
CAGGTTA  (ornithine  decarboxylase);  and  5'-CCAGAGCAAGAGAGGT 
ATCC  and  5'-CTGTGGTGGTGAAGCTGTAG  ((3-actin).  RT-PCRs  were 
performed  with  One  Step  RT-PCR  kits  (Invitrogen),  and  the  number  of  cycles 
for  each  reaction  template  was  optimized  to  ensure  that  all  reactions  being 
compared  were  within  the  linear  portion  of  the  amplification  curve.  Reaction 
products  were  analyzed  on  1.25%  agarose  gels,  stained  with  ethidium  bromide, 
and  visualized  under  UV  light.  Control  reactions  were  done  without  reverse 
transcriptase  (with  the  thermostable  DNA  polymerase  alone)  to  ensure  that 
there  was  no  contaminating  DNA  in  the  RNA  samples  being  assayed. 

Cell  death  assays.  Measurement  of  necrosis  and  apoptosis  within  RAW  264.7 
cells  was  done  with  the  Cell  Death  Detection  ELISA  Plus  (Roche)  according  to 
the  manufacturer’s  instructions.  Briefly,  104  cells  per  sample  were  cultured  in  a 
96-well  plate,  and  B.  anthracis  34F2  spores,  camptotheein,  d,i  -a-difluoromelh- 
ylornithine  (DFMO),  or  a  combination  of  these  were  added  in  fresh  medium. 
Untreated  cells  were  also  given  fresh  medium  without  any  added  constituents.  At 
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the  appropriate  lime  points,  the  presence  of  histone-associated  DNA  fragments 
was  assayed  in  both  the  culture  supernatant  and  the  cytoplasm  of  infected  cells. 
As  prescribed  by  the  manufacturer,  we  assigned  a  value  of  1.0  to  the  levels  of 
histone-associated  DNA  detected  in  untreated  control  cells,  and  the  relative 
values  of  experimental  samples  were  recorded  as  a  ratio  to  the  control.  Exper¬ 
iments  were  performed  at  least  three  times  each,  and  within  each  experiment 
measurements  were  taken  from  at  least  three  replicate  samples.  Statistical  sig¬ 
nificance  was  measured  by  Student's  two-tailed  /  test.  Although  different  exper¬ 
iments  varied  as  to  the  overall  magnitude  of  cell  death  recorded  (perhaps  re¬ 
flecting  variations  in  the  age  of  RAW  264.7  cells),  the  results  depicted  in  Fig.  3 
are  typical,  and  in  every  instance  the  differences  between  untreated  cells  and 
infected  cells,  camptothecin-treated  cells,  and  infected  cells  treated  with  DFMO 
were  statistically  significant  (P  <  0.05). 

Thiazolyl  blue  assays  for  measurement  of  cell  viability.  Killing  of  RAW264.7 
cells  by  B.  anthracis  lethal  toxin  was  assayed  with  the  thiazolyl  blue  (Sigma 
M5655)  assay.  Before  assay,  the  medium  was  changed  to  DMEM  plus  10  mM 
HOEPES,  pH  7.4  (DMEM-HEPES)  and  the  cell  concentration  adjusted  to  4  x  105 
cells/ml.  Cells  (6  x  104  per  well)  were  plated  in  a  96-well  plate  and  allowed  to 
settle  at  room  temperature  for  30  min.  Cells  were  incubated  at  37°C  and  5%  C02 
for  16  h.  Medium  was  aspirated  and  replaced  with  DMEM-HEPES  containing 
test  toxin  concentrations  in  triplicate.  Lethal  toxin-containing  samples  ranged 
from  4  to  500  ng/ml  of  lethal  factor  (recombinant  lethal  factor  and  protective 
antigen;  List  Biological  Laboratories)  plus  a  constant  concentration  (0.5  jxg/ml) 
of  protective  antigen,  with  or  without  a  constant  concentration  (0.5  p.g/ml)  of 
lipopolysaccharide  (LPS;  from  Salmonella  enieritidis.  Sigma  L-4774).  Plates  were 
incubated  for  4  h  at  37<5C,  5%  C02,  and  50  p,l  of  thiazolyl  blue  dissolved  in 
DMEM-HEPES  at  1.5  mg/ml  was  added  to  each  well  and  incubated  for  L  h  at 
37°C  and  5%  C02.  Aspirated  medium  was  then  replaced  with  100  gU  of  0.5% 
sodium  dodecyl  sulfate-25  mM  HCI  in  90%  isopropanol  and  vortexed  gently,  and 
absorption  was  read  at  570  nm  to  650  nm  on  a  Molecular  Devices  SpectraMAX 
250  microplate  reader.  Lysis  of  RAW  264.7  cells  was  also  measured  with  the 
same  range  of  lethal  toxin  and  lethal  toxin  plus  LPS  with  the  trypan  blue  visual 
assay.  Briefly,  cells  were  grown  and  treated  exactly  as  above,  but  scraped  after  4  h 
and  mixed  in  a  1:1  mixture  of  0.1%  trypan  blue  to  cell  suspension,  counted  on  a 
hemacytometer,  and  percent  lysis  was  calculated  by  percentage  of  blue  cells. 

Data  availability.  The  expression  data  collected  in  this  study  are  publicly 
available  via  theTIGR  website  (http://pga.tigr.org)  and  the  ArrayExpress  data¬ 
base  (accession  number  E-TIGR-84).  Complete  lists  of  genes  that  were  identi¬ 
fied  as  differentially  expressed  in  response  to  any  of  the  experimental  conditions 
described  here  are  in  the  supplemental  material  at  http://www.sitemaker.umich 
.edu/hanna  lab 


RESULTS  AND  DISCUSSION 

Although  many  of  the  individual  virulence  factors  employed 
by  B .  anthracis  in  causing  disease  are  known,  the  effects  of  the 
infecting  bacteria  on  host  gene  expression  have  yet  to  be 
clearly  defined.  It  was  therefore  of  interest  to  examine  the 
transcriptional  changes  that  occur  within  infected  macro¬ 
phages  and  identify  genes  that  are  induced  or  repressed  in 
response  to  B.  anthracis.  We  also  sought  to  compare  these 
trends  to  those  observed  in  macrophages  exposed  to  lethal 
toxin,  as  well  as  to  cells  that  had  been  exposed  to  Salmonella 
lipopolysaccharide  (LPS),  a  reagent  that  is  known  to  induce  a 
widespread  activation  program  within  macrophages  (46). 

The  RAW  264.7  murine  macrophage-like  cell  line  provided 
a  convenient  cell  culture  model  system  for  these  experiments, 
as  it  has  been  used  widely  in  the  past  for  studies  of  B.  anthracis 
as  well  as  other  bacterial  pathogens  (12,  23,  26,  30,  56).  Initial 
experiments  matched  previously  published  work  (12,  48)  and 
showed  that  1  h  after  infection  vegetative  bacilli  could  be  seen 
inside  or  closely  associated  with  the  macrophages  (Fig.  1,  panel 
B).  After  3  h  the  bacteria  were  beginning  to  multiply  (panel  C), 
and  6  h  after  infection  the  majority  of  macrophages  had  lysed, 
and  bacteria  were  growing  in  the  extracellular  medium  (panel 
D). 


Microarray  analysis  of  macrophages  infected  with  B.  anlhra - 
cis  endospores.  In  order  to  define  the  patterns  of  gene  expres¬ 
sion  associated  with  these  events,  we  isolated  total  RNA  from 
RAW  264.7  cells  following  infection  with  B.  anthracis  Sterne 
34F2  spores.  Samples  were  harvested  every  hour  for  6  h,  after 
which  relatively  few  intact  murine  cells  were  observed,  and 
RNA  yields  decreased  dramatically.  RNA  samples  were  used 
to  make  cDNA  probe  mixtures,  which  were  then  applied  to 
DNA  microarrays  containing  target  sequences  corresponding 
to  27,648  murine  genes  and  expressed  sequence  tags.  For  com¬ 
parison,  mRNA  samples  from  parallel  mock  infections,  in 
which  B.  anthracis  spores  were  omitted  and  only  fresh  medium 
was  added  to  the  cells,  were  also  isolated  and  prepared  in  the 
same  way.  In  each  case  data  were  collected  and  processed  as 
described  in  Materials  and  Methods,  and  genes  that  were  dif¬ 
ferentially  expressed  in  a  given  sample  were  selected  at  a  99% 
confidence  level. 

Prior  to  analyzing  gene  expression  in  infected  mouse  mac¬ 
rophages,  it  was  necessary  to  determine  to  what  extent  gene 
expression  changed  over  time  in  the  mock-infected  cells,  as 
these  data  provide  an  indication  of  the  random  fluctuation  of 
gene  expression  levels  within  our  experiment.  Accordingly,  we 
compared  the  expression  profiles  of  cells  at  0,  3,  and  6  h  after 
a  mock  infection  in  which  the  cells  were  given  fresh  medium 
(as  in  all  other  experiments)  without  challenge.  Within  this 
time  interval,  gene  expression  in  RAW  264.7  macrophages 
changed  very  little:  between  0  and  3  h  after  the  mock  infection, 
the  median  change  in  expression  level  was  1.0035-fold,  with 
95%  changing  by  less  than  1.26-fold  and  only  six  genes  chang¬ 
ing  more  than  2-fold  (supplemental  Fig.  SI).  Between  0  and 
6  h,  the  overall  level  of  change  was  even  lower,  with  median 
change  of  1.0025-fold,  95%  changing  by  less  than  1. 21-fold, 
and  only  two  genes  changing  in  expression  level  by  twofold  or 
more.  These  measurements  imply  that  changes  observed  in 
comparing  different  treatments,  or  cells  treated  the  same  way 
at  different  times,  are  likely  to  be  due  to  the  experimental 
difference  and  not  due  to  stochastic  changes  in  gene  expres¬ 
sion. 

Previous  studies  with  this  experimental  system  and  our  own 
results  have  shown  that  roughly  3  h  after  inoculation,  B.  an¬ 
thracis  is  able  to  escape  from  the  phagolysosome  into  the  host 
cell  cytoplasm.  This  event  thus  separates  the  early  phase  of 
infection  (those  events  occurring  in  the  phagolysosome)  from 
the  late  phase  (the  events  occurring  within  the  host  cell  cyto¬ 
plasm,  leading  up  to  killing  of  the  phagocyte  and  escape  into 
the  extracellular  environment). 

In  our  examination  of  macrophage  gene  expression  during 
infection  by  B.  anthracis ,  we  analyzed  these  stages  separately 
and  began  by  assessing  the  gene  expression  patterns  that  occur 
early  during  infection,  which  encompasses  phagocytosis  of  en¬ 
dospores,  germination,  and  the  beginning  of  bacterial  replica¬ 
tion.  Within  this  time  frame,  493  genes  met  our  criteria  for 
significant  differential  expression  (see  Materials  and  Methods 
for  descriptions  of  these  criteria),  with  213  induced  and  280 
repressed  in  at  least  one  of  the  three  time  points  (see  supple¬ 
mental  Table  SI).  These  genes  were  sorted  by  functional  fam¬ 
ily  with  the  Onlo-Express  software  package,  which  uses  the 
gene  ontology  annotation  system  to  sort  genes  by  pathway  and 
function  (14).  This  allowed  the  identification  of  biological  pro¬ 
cesses  that  were  related  to  the  genes  within  a  given  data  set, 
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FIG.  1.  Infection  of  RAW  264.7  cells  by  Bacillus  anthracis.  Murine  macrophage-like  RAW  264.7  cells  are  shown  untreated  (A)  and  1  h  (B),  3  h 
(C),  and  6  h  (D)  after  inoculation  with  B.  anthracis  34F2  spores  at  a  multiplicity  of  infection  of  10:1. 


and  an  assessment  of  the  statistical  significance  of  any  over¬ 
representations  that  were  found.  Those  biological  processes 
with  a  statistically  significant  overrepresentation  (P  <  0.05) 
and  at  least  three  members  present  in  the  early  (1  to  3  h 
postinoculation)  data  sets  are  shown  in  Table  1. 

Gene  expression  during  early  infection  by  B.  anthracis  was 
marked  by  several  distinct  trends.  Perhaps  the  most  prominent 
of  these  is  the  induction  of  both  innate  immune  response  and 
stress-related  genes.  These  loci  include  the  CD14  gene,  which 
encodes  the  receptor  for  the  complex  of  bacterial  lipopolysac- 
charide  (LPS)  and  LPS-binding  protein,  and  has  been  recently 
and  previously  shown  to  be  essential  for  the  innate  immune 
response  to  bacterial  infection  (53).  The  product  of  this  gene 
has  also  been  recently  shown  to  mediate  the  innate  immune 
response  to  gram-positive  bacteria  by  binding  to  cell  wall- 
derived  peptidoglycan-polysaccharide  complexes,  and  this  may 
explain  its  induction  in  response  to  B.  anthracis  (35).  Also 


observed  was  an  induction  of  genes  associated  with  Toll-like 
receptor  signaling  (e.g.,  I-kB  zeta,  which  has  been  seen  in  other 
studies  to  be  rapidly  induced  in  response  to  activation  of  the 
Myd88-dependent  signaling  pathway,  commonly  found  down¬ 
stream  of  different  Toll-like  receptors)  (19,  34)  as  well  as  loci 
previously  known  to  play  a  role  in  p38  mitogen-activated  pro¬ 
tein  kinase  activation  (e.g.,  Mkk3,  Mkklipl,  phospholipase  A2, 
and  c-Myc),  which  is  the  pathway  best  characterized  as  the 
target  for  disruption  by  anthrax  lethal  toxin  (7,  41). 

It  seems  clear  from  the  functional  families  represented  in 
both  the  induced  and  the  repressed  data  sets  that  the  first  3  h 
following  infection  by  B.  anthracis  spores  are  marked  by  major 
shifts  in  cellular  signaling,  as  evidenced  by  the  changes  in  the 
expression  levels  of  genes  involved  in  receptor-linked  signal 
transduction,  ion  transport,  protein  phosphorylation  and  de¬ 
phosphorylation,  and  GTP-linked  signaling  cascades.  There 
were  also  significant  changes  in  genes  previously  linked  to 
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TABLE  1.  Functional  families  overrepresented  in  the  macrophage  transcriptional  response  to  early  stages  of  B.  anthracis  infection 


Induced 

Repressed 

Biological  process 

P 

Biological  process 

P 

Biogenic  amine  metabolism 

LI  E-4 

Protein  biosynthesis 

0.0351 

Innate  immune  response  (inflammatory 

1.2  E-4 

Chromosome  organization  and  biogenesis 

0.0429 

response) 

Receptor-linked  signal  transduction 

1.3  E-4 

Response  to  stress  (pathogen/DNA 

1.8  E-4 

Protein  folding 

6.5  E-4 

damage/wounding) 

Protein  biosynthesis 

0.0038 

Cell  migration 

2.2  E-4 

Cell  homeostasis 

0.0053 

Regulation  of  apoptosis 

7.7  E-4 

Regulation  of  transcription 

0.0092 

Cytoskelelon  organization  and  biogenesis 

0.0013 

Regulation  of  cel!  cycle 

0.0100 

Sterol  metabolism 

0.0028 

Ion  transport 

0.0145 

Cell  growth  and/or  maintenance 

0.0048 

Protein  kinase  signaling  cascades 

0.0158 

DNA  repair 

0.0088 

Proteolysis 

0.0253 

Regulation  of  cell  cycle 

0.0151 

Defense  response  to  pathogenic  bacteria 

0.0347 

Receptor-linked  signal  transduction 

0.0164 

Protein  phosphorylation  and 

0.0395 

Amino  acid  metabolism 

0.0183 

dephosphorylation 

Cell-cell  adhesion 

0.0249 

Regulation  of  apoptosis 

0.0439 

GTPase-mediated  signal  transduction 
Morphogenesis 

0.0264 

0.0265 

Cell-cell  adhesion 

0.0443 

control  of  cell  cycle  (e.g.,  c-Myc),  regulation  of  apoptosis  (e.g., 
Bel- 10),  and  cellular  growth  and  homeostasis  (e.g.,  Rabl).  Col¬ 
lectively,  these  transcriptional  patterns  likely  reflect  the  acti¬ 
vation  of  the  murine  macrophages  following  phagocytosis  of 
the  B.  anthracis  spores.  Concurrent  with  these  changes,  we 
noted  expression  level  changes  in  loci  involved  in  cell  migra¬ 
tion  (e.g.,  somatostatin),  cytoskeletal  biogenesis  and  organiza¬ 
tion  (e.g.,  ^-tubulin),  and  cell  adhesion  (e.g.,  protocadherin 
10),  which  may  reflect  the  fact  that  in  each  of  the  three  forms 
of  anthrax,  phagocytosis  and  activation  of  the  macrophage  are 
closely  followed  by  migration  of  the  phagocyte  to  the  regional 
lymph  nodes. 

The  lymph  nodes  are  the  site  where  the  late  stages  of  B. 
anthracis  infection  begin,  and  in  cell  culture  this  time  frame  is 
marked  by  the  escape  of  the  bacteria  from  the  phagolysosome 
into  the  host  cell  cytoplasm,  replication  within  the  macro¬ 
phage,  and  eventual  killing  of  the  host  cell  with  subsequent 
release  of  the  bacterium  into  the  extracellular  environment. 
These  events  were  observed  in  previous  work  and  in  our  own 
experiments  to  occur  roughly  4  to  6  h  postinoculation  (Fig.  1), 
In  order  to  identity  the  transcriptional  patterns  underlying  the 
macrophage  response  to  the  later  stages  of  infection,  we  ana¬ 
lyzed  total  RNA  samples  taken  4,  5,  and  6  h  postinoculation  in 
the  same  way  the  studies  of  the  early  stages  of  infection  had 
been  done.  Within  this  time  frame,  480  genes  showed  signifi¬ 
cant  differential  expression,  with  292  genes  induced  and  188 
genes  repressed  in  at  least  one  of  the  three  time  points  exam¬ 
ined  (see  supplemental  Table  S2). 

With  the  methods  described  above,  differentially  expressed 
genes  were  sorted  according  to  functional  category,  and  bio¬ 
logical  processes  whose  components  showed  statistically  signif¬ 
icant  (P  <  0.05)  overrepresentation  are  listed  in  Table  2.  Sa¬ 
lient  findings  included  increases  in  the  mRNA  levels  of  many 
genes  associated  with  the  immune  response  and  cytokine  sig¬ 
naling  (e.g.,  macrophage  colony-stimulating  factor,  tumor  ne¬ 
crosis  factor  alpha,  and  GADD45-(3),  reflecting  the  growing 
severity  of  the  infection.  Many  of  the  immune  response  genes 
that  were  upregulated  in  early  infection  continued  to  increase 
in  expression  level  into  the  later  stages  of  the  infection,  and  a 


variety  of  genes  generally  linked  to  the  immune  response  were 
induced  for  the  first  time  in  this  time  frame.  In  general,  the 
immune  response  genes  differentially  expressed  during  late 
infection  by  B.  anthracis  infection  were  the  same  genes  iden¬ 
tified  as  part  of  the  stereotypical  innate  immune  response 
transcriptional  program.  This  set  of  transcriptional  responses  is 
triggered  by  a  range  of  bacterial  pathogens  and  has  recently 
been  characterized  by  several  groups  (4,  6,  38).  One  notable 
parallel  is  in  the  induction  of  NF-kB  and  a  range  of  associated 
genes  and  regulatory  targets  (e.g.,  I-kB,  I-kB  kinase)  which 
have  been  shown  in  several  other  studies  to  play  pivotal  roles 
in  immune  system  activation  (3,  4,  6,  38). 

In  addition  to  the  “general  immune  response  genes,”  we 
continued  to  observe  differential  regulation  of  genes  involved 
in  regulation  of  cell  cycle,  intra-  and  extracellular  signaling, 
and  adhesion  both  to  other  cells  and  to  the  extracellular  ma¬ 
trix.  Many  of  these  genes  were  identified  early  in  infection,  and 
it  would  appear  that  many  of  the  processes  begun  soon  after 
phagocytosis  are  still  occurring  late  in  B.  anthracis  infection. 
Intriguingly,  the  most  significantly  overrepresented  functional 
family  was  the  set  of  genes  involved  in  regulating  apoptosis. 
Several  of  these  genes  were  observed  to  be  differentially  ex¬ 
pressed  during  early  infection,  but  the  trend  was  much  more 
pronounced  4  to  6  h  after  inoculation,  and  this  particular 
functional  family  included  several  of  the  most  highly  induced 
genes  identified  (e.g.,  ornithine  decarboxylase  and  the  orphan 
nuclear  receptor  NURR1,  both  upregulated  ~  14-fold  late  in 
infection).  Significantly,  the  time  frame  in  which  the  apoptosis- 
associated  genes  became  most  highly  induced  was  essentially 
the  same  in  which  we  first  observed  cell  death  occurring  in 
infected  macrophages  (Fig.  ID),  suggesting  that  the  cell  death 
caused  by  intracellular  B.  anthracis  is  apoptotic  in  nature.  This 
is  consistent  with  a  recent  study  by  Popov  et  al.  (45),  which 
reported  membrane  permeability  changes  commonly  associ¬ 
ated  with  early  stages  of  apoptosis  in  B.  anthracis- infected 
macrophages,  as  well  as  with  several  other  studies  that  have 
shown  that  anthrax  lethal  toxin,  a  binary  exotoxin  which  is 
required  for  full  virulence  in  B.  anthracis ,  is  able  to  efficiently 
induce  apoptosis  and/or  necrosis  in  a  variety  of  cell  types  by 
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TABLE  2.  Functional  families  overrepresented  in  the  macrophage  transcriptional  response  to  late  stages  of  B.  anthracis  infection 


Induced 

Repressed 

Biological  process 

P 

Biological  process 

P 

Regulation  of  apoptosis 

2.6  E-6 

Cell-matrix  adhesion 

5.0  E-5 

Regulation  of  transcription 

1.0  E-5 

G-protein  coupled  receptor  protein 

3.1  E-4 

signaling  pathway 

Nitrogen  metabolism 

2.0  E-5 

Angiogenesis 

0.0013 

Immune  response  (inflammatory 

2.0  E-5 

Regulation  of  cell  cycle 

0.0100 

response) 

Regulation  of  cell  cycle 

5.0  E-5 

Cell  adhesion 

0.0109 

Intracellular  signaling  cascade 

8.0  E-5 

Regulation  of  apoptosis 

0.0113 

Protein  phosphorylation  and 

9.0  E-5 

Immune  response 

0.0175 

desphosphorylation 

Cytokine  biosynthesis 

1.2  E-4 

Intracellular  signaling  cascade 

0.0212 

Amine  catabolism 

3.8  E-4 

Cytoskeleton  organization  and 

0.0012 

biogenesis 

Protein  kinase  cascade 

0.0025 

Cell  differentiation 

0.0026 

Cell-matrix  adhesion 

0.0034 

Cell  adhesion 

0.0039 

Receptor-linked  signal  transduction 

0.0094 

GTPase-mediated  signaling  cascade 

0.0182 

Proteolysis 

0.0345 

inactivating  a  p38  mitogen -activated  protein  kinase  signaling 
pathway  (15,  33,  41,  44,  45). 

Confirmation  of  specific  gene  expression  trends  by  RT-PCR. 

In  general,  our  microarray  analysis  of  B.  anthracis  infection 
suggested  that  many  of  the  transcriptional  patterns  character¬ 
istic  of  a  general  macrophage  response  to  bacterial  pathogens 
are  also  found  in  the  response  to  B.  anthracis.  In  order  to  more 
completely  assess  the  reliability  of  the  arrays  and  the  data 
derived  from  them,  we  sought  to  validate  the  RNA  expression 
trends  by  an  alternative  method.  To  do  this,  we  selected  rep¬ 
resentative  genes  from  each  of  three  broad  categories  (induced 
following  expression;  repressed  following  infection;  and  un¬ 
changed  during  infection)  and  repeated  our  measurements  of 
RNA  expression  with  RT-PCR  (Fig.  2).  Total  RNA  samples 
were  isolated  as  in  the  case  of  microarray  analysis,  and  RT- 
PCRs  were  performed  with  equivalent  amounts  of  starting 
RNA.  Preliminary  experiments  were  done  to  determine  the 
optimal  number  of  cycles  for  each  PCR  template  such  that  all 
of  the  reactions  being  compared  were  within  the  linear  portion 
of  the  amplification  curve.  The  increase  in  the  level  of  orni¬ 
thine  decarboxylase  expression  that  we  observed  by  microarray 
following  infection  was  confirmed  in  our  RT-PCR  analysis,  as 
was  the  decrease  in  the  level  of  methylmalonate  semialdehyde 
dehydrogenase  expression.  The  gene  encoding  p-actin,  which 
did  not  show  an  appreciable  change  in  expression  level  in  our 
microarray  measurements  and  thus  was  an  appropriate  con¬ 
trol,  similarly  remained  at  a  relatively  constant  level  of  expres¬ 
sion  throughout  the  experimental  time  course  when  measured 
by  RT-PCR. 

Confirmation  of  B.  anthracis- induced  apoptosis  and  identi¬ 
fication  of  a  major  antiapoptotic  factor  induced  in  infected 

macrophages.  Intuitively,  it  seems  that  the  ability  to  cause 
apoptosis  in  host  macrophages  would  almost  certainly  be  an 
advantage  for  B.  anthracis ,  as  apoptotic  cells  are  much  less 
likely  than  necrotic  cells  to  burst  and  release  cytokines  and 
other  proinflammatory  signals  into  the  extracellular  environ¬ 


ment.  Since  our  data  as  well  as  those  of  other  groups  indicate 
that  B.  anthracis  activates  the  host  macrophage  during  infec¬ 
tion,  it  may  be  that  induction  of  apoptosis  is  a  means  by  which 
the  bacterium  is  able  to  kill  macrophages  silently  and  to  a  large 
extent  delay  or  prevent  a  heightened  immune  response  in  the 
area  of  infection. 

Because  the  microarray  data  seemed  to  point  strongly  to¬ 
ward  apoptosis  as  the  underlying  mode  of  intracellular  B.  an- 
thracis- induced  cell  death,  we  sought  to  confirm  this  possibility 
through  other  means.  Accordingly,  we  infected  RAW  264.7 
cells  with  B.  anthracis  Sterne  34F2  spores  and  used  a  commer¬ 
cially  available  enzyme-linked  immunosorbent  assay  kit  (Cell 
Death  Detection  kit;  Roche)  to  detect  and  differentiate  be¬ 
tween  apoptosis  and  necrosis.  Specifically,  the  kit  detects  his¬ 
tone-associated  DNA  fragments,  which  remain  membrane 
bound  following  apoptosis  but  are  released  into  the  extracel¬ 
lular  environment  during  necrosis.  The  two  forms  can  thus  be 
measured  independently,  and  the  precise  mode  of  cell  death 
that  is  occurring  can  be  determined. 

Microscopic  observation  of  infected  macrophages  showed 
that  3  h  after  inoculation  with  B.  anthracis  Sterne  34F2  spores, 
bacteria  were  replicating  within  the  host  cells  but  had  not 
escaped,  while  after  6  h  there  was  appreciable  cell  death  and 
bacteria  could  be  seen  in  the  extracellular  medium  (Fig.  1, 
panels  C  and  D).  Our  cell  death  assays  closely  mirrored  these 
observations  and  detected  essentially  no  cell  death  3  h  posti¬ 
noculation.  We  did,  however,  detect  significant  apoptosis  6  h 
after  the  beginning  of  infection  (Fig.  3).  The  levels  of  apoptosis 
in  B.  anthracis- infected  cells  at  the  final  time  point  were  similar 
to  the  levels  detected  in  cells  treated  with  camptothecin,  a 
well-characterized  inducer  of  apoptosis  (27). 

As  noted  earlier,  many  of  the  most  highly  induced  genes  late 
in  infection  are  related  to  the  regulation  of  apoptosis.  One 
particularly  interesting  example  is  the  locus  encoding  the  en¬ 
zyme  ornithine  decarboxylase,  which  was  induced  by  approxi¬ 
mately  15-fold  late  in  the  infection.  This  enzyme  catalyzes  a 
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FIG.  2.  RT-PCR  measurement  of  expression  levels  for  representative  genes.  A  Expression  levels  for  selected  genes  as  measured  by  DNA 
microarray.  Levels  shown  reflect  induction  or  repression  in  cells  3  and  6  h  postinfection  relative  to  untreated  cells.  B  RT-PCR  analysis  of  gene 
expression.  RT-PCRs  were  performed  on  equal  amounts  of  starting  RNA,  and  control  reactions  (not  shown)  demonstrated  that  there  was  no 
contaminating  DNA  present  in  the  RNA  samples. 


key  step  in  the  synthesis  of  spermine  and  spermidine,  key 
polyamines  that  are  known  to  be  essential  for  cell  proliferation 
and  differentiation.  Significantly,  these  molecules  also  exert  a 
powerful  protective  effect  against  cell  death  induced  by  a  va¬ 
riety  of  cellular  stresses  (8,  32),  and  a  recent  study  by  Park  et 
al.  (40)  showed  that  overexpression  of  ornithine  decarboxylase 
is  sufficient  to  protect  cells  from  apoptosis  induced  by  H202, 
7-irradiation,  and  a  variety  of  chemotherapeutic  agents.  Since 
it  appeared  that  induction  of  apoptosis  may  be  a  key  element 
in  B.  anthracis  pathogenesis,  and  elevated  levels  of  ornithine 
decarboxylase  expression  have  not  been  reported  in  studies 
examining  the  macrophage  response  to  a  variety  of  other  bac¬ 
terial  pathogens  (3,  4,  6,  38),  we  sought  to  examine  the  role  of 
ornithine  decarboxylase  in  more  detail  and  determine  the  ex¬ 
tent  of  its  importance  in  infection. 


Ornithine  decarboxylase  expression  has  long  been  known  to 
be  a  marker  associated  with  a  variety  of  cancers,  and  several 
chemotherapeutic  agents  are  able  to  specifically  inhibit  orni¬ 
thine  decarboxylase.  The  most  specific  and  best  characterized 
of  these  is  n,L-a-difluoromethylornithine  (DFMO),  which  is  an 
irreversible  suicide  inhibitor  of  ornithine  decarboxylase.  Con¬ 
veniently,  this  compound  has  been  used  in  previous  studies 
examining  ornithine  decarboxylase  activity  in  RAW  264.7  cells 
and  has  been  shown  to  effectively  inhibit  the  enzyme  without 
affecting  cell  viability  (21).  We  therefore  took  advantage  of 
these  properties  and  used  DFMO  to  inhibit  ornithine  decar¬ 
boxylase  in  B.  anthracis- infected  RAW  cells  in  order  to  better 
understand  the  role  of  this  enzyme  in  the  response  to  infection. 

Using  the  enzyme-linked  immunosorbent  assay  system  to 
assess  cell  death,  we  measured  necrosis  and  apoptosis  in  mac- 
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FIG.  3.  Cell  death  (necrosis  and  apoptosis)  in  B.  anthracis  infection  as  measured  by  enzyme-linked  immunosorbent  assay.  Measurements  of 
necrosis  and  apoptosis  have  been  normalized  in  each  case  to  measurements  obtained  from  untreated  cells,  such  that  levels  of  cell  death  are 
expressed  as  the  ratio  of  experimental  to  untreated  cells.  Data  shown  are  from  representative  experiments,  and  each  column  is  the  average  of  at 
least  three  individual  measurements.  Error  bars  indicate  the  standard  deviation  of  each  data  set. 


ro phages  that  were  infected  with  B.  anthracis  spores,  treated 
with  DFMO,  or  infected  in  the  presence  of  DFMO,  and  com¬ 
pared  them  to  untreated  cells  and  cells  treated  with  campto- 
thecin  (Fig.  3).  Cells  treated  with  inhibitor  alone  were  essen¬ 
tially  indistinguishable  from  untreated  cells,  indicating  that,  as 
previous  studies  had  shown,  DFMO  did  not  have  a  discernible 
effect  on  RAW  264.7  cells.  In  comparing  infected  cells  with 
those  that  had  been  infected  in  the  presence  of  DFMO,  we 
noted  a  significant  difference  in  how  rapidly  apoptosis  was 
induced.  Infected  cells  treated  with  DFMO  showed  signs  of 
apoptosis  3  h  after  infection  began,  while  at  this  time  point 
cells  infected  without  the  inhibitor  did  not  show  significant  cell 
death.  These  data  were  consistent  with  microscopic  observa¬ 
tions,  which  showed  infection  progressing  more  rapidly  in  cells 
treated  with  inhibitor  (data  not  shown),  and  suggest  that  orni¬ 
thine  decarboxylase  is  a  key  suppressor  of  apoptosis  during 
infection  by  B.  anthracis.  Together,  these  findings  further  high¬ 
light  the  importance  of  apoptosis  induction  as  an  element  of  B. 
anthracis  pathogenesis  and  are  consistent  with  a  recent  study 
that  suggested  that  small-molecule  inhibitors  of  apoptosis  may 
be  therapeutically  useful  in  treating  anthrax  and  in  preventing 
the  disease  in  those  exposed  to  spores  (43). 

Finally,  it  is  of  interest  to  point  out  that  ornithine  decarbox¬ 
ylase  expression  has  been  examined  in  response  to  a  variety  of 
bacterial  pathogens,  including  enterohemorrhagic  E.  coli,  Sal¬ 
monella  enterica  serovar  Typhi,  Listeria  monocytogenes ,  Myco¬ 
bacterium  tuberculosis ,  Bordetella  pertussis ,  and  Staphylococcus 
aureus ,  and  in  each  case  no  discernible  induction  of  expression 
was  observed.  It  is  tempting  to  speculate  that  a  detectable 
increase  in  expression  of  this  gene  might  be  unique  to  B. 
anthracis  and  thus  may  serve  as  a  useful  marker  for  specific 
diagnosis  of  early-stage  anthrax,  though  we  note  that  the  mac¬ 
rophage  responses  to  close  relatives  like  Bacillus  cereus  as  well 
as  other  to  biowarfare  agents  have  not  yet  been  characterized, 
and  these  data  will  be  necessary  to  fully  test  this  hypothesis. 

Gene  expression  patterns  induced  by  anthrax  lethal  toxin  in 
macrophages.  Because  anthrax  lethal  toxin  is  known  to  have 


profound  effects  on  host  macrophages,  including  disruption  of 
a  variety  of  signaling  pathways  and  induction  of  apoptosis  (20, 
25),  it  was  of  interest  to  determine  to  what  extent  the  gene 
expression  patterns  observed  in  infected  cells  could  be  attrib¬ 
uted  to  the  actions  of  lethal  toxin  alone.  This  toxin  consists  of 
the  two  proteins  protective  antigen,  which  facilitates  entry  into 
host  cells,  and  lethal  factor,  which  is  a  zinc  metalloprotease 
that  cleaves  members  of  the  MEK  signaling  family.  It  is  known 
to  disrupt  a  variety  of  signaling  pathways  within  the  host  cell, 
and  previous  studies  have  shown  that  the  effects  of  lethal  toxin 
differ  depending  on  the  activation  state  of  the  macrophage  (37, 
41).  We  therefore  examined  gene  expression  in  both  activated 
and  unactivated  macrophages  following  treatment  with  an¬ 
thrax  lethal  toxin  and  compared  these  expression  patterns  to 
those  observed  in  untreated  cells  as  well  as  those  that  were 
activated  but  not  treated  with  toxin. 

RAW  264.7  cells  were  treated  with  lethal  toxin,  Salmonella 
lipopolysaceharide  (LPS),  or  a  combination  of  the  two,  and 
total  RNA  was  prepared  as  described  for  infected  cells.  LPS 
was  used  at  a  constant  concentration  of  500  ng/ml,  and  at  this 
level  it  has  been  previously  shown  to  induce  a  robust,  general 
activation  program  within  RAW  264.7  cells  without  signifi¬ 
cantly  affecting  cell  viability  (46).  Lethal  toxin  was  used  at  a 
concentration  (8  ng/ml)  that  previous  studies  and  our  own 
preliminary  experiments  had  shown  was  slightly  below  the 
threshold  at  which  cell  death  began  to  occur  (data  not  shown). 
RNA  was  collected  3  and  6  h  after  addition  of  toxin  and/or 
LPS,  and  samples  were  analyzed  and  differentially  expressed 
genes  were  selected  at  the  99%  confidence  level  as  described 
for  infection. 

As  expected,  the  expression  patterns  observed  in  cells 
treated  with  Salmonella  LPS  were  reflective  of  a  profound 
change  in  cellular  signaling,  and  in  general  these  responses 
were  somewhat  more  intense  than  those  observed  in  any  other 
experimental  condition,  in  that  the  changes  in  both  individual 
genes  (see  supplemental  Table  S3)  and  functional  families 
(Table  3)  showed  higher  levels  of  statistical  significance.  The 
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TABLE  3.  Functional  families  overrepresented  in  the  macrophage  transcriptional  response  to  treatment  with  Salmonella 

lipopolysaccharide  (500  ng/ml) 


Induced 

Repressed 

Biological  process 

P 

Biological  process 

P 

Cell  growth 

8.9  E-8 

Mitotic  cell  cycle.  S  phase 

4.5  E-8 

Regulation  of  cell  cycle 

5.2  E-7 

Regulation  of  cell  cycle 

4.3  E-7 

Lymph  gland  development 

1.0  E-5 

Cytokinesis 

1.3  E-5 

Regulation  of  apoptosis 

2.0  E-5 

Lipid  metabolism 

8.0  E-5 

DNA  repair 

9.2  E-5 

Cell  motility 

3.2  E-4 

Immune  response 

1.9  E-4 

Protein  phosphorylation 

2.1  E-4 

Angiogenesis 

2.7  E-4 

Mitotic  cell  cycle,  M  phase 

7.5  E-4 

Regulation  of  transcription 

4.5  E-4 

Endocytosis 

0.0032 

Intracellular  signaling  cascade 

2.6  E-4 

Apoptosis 

0.0233 

Cell  differentiation 

0.0021 

Protein  folding 

0.0074 

Amine  metabolism 

0.0026 

Regulation  of  transcription 

0.0031 

Cell  proliferation 

0.0055 

Intracellular  signaling 

0.0066 

Stress  response  (to  pathogen/ 

DNA  damage/wounding) 

0.0130 

cascade 

genes  identified  in  this  set  correspond  well  to  those  identified 
in  previous  studies  examining  expression  trends  in  macro¬ 
phages  treated  with  a  variety  of  bacterial  pathogens  (4,  6,  38, 
46)  and  include  many  classical  proteins  involved  in  immune 
signaling  (e.g.,  tumor  necrosis  factor  alpha,  interleukin-la, 
macrophage  colony-stimulating  factor,  and  NF-kB). 

With  the  expression  patterns  from  both  untreated  cells  and 
LPS-activated  cells  as  control  baselines,  we  examined  gene 
expression  in  both  groups  after  treatment  with  lethal  toxin. 
Surprisingly,  the  cells  treated  with  lethal  toxin  alone  (unacti¬ 
vated  macrophages)  showed  essentially  no  significant  changes 
in  gene  expression  compared  to  untreated  cells.  The  expres¬ 
sion  profiles  from  these  cells  and  those  from  untreated  control 
cells  (described  above)  looked  very  similar,  with  only  17  genes 
changing  by  ^2-fold  throughout  the  6-h  time  course  and  no 
coherent  pattern  with  regard  to  function  within  this  small 
group  (supplemental  Table  S4).  The  finding  that  gene  expres¬ 
sion  did  not  change  appreciably  throughout  a  toxin  treatment 
that  caused  appreciable  cell  death  (10  to  20%)  was  somewhat 
surprising,  but  it  matched  the  observations  in  a  similar  study 
recently  published  by  Tucker  et  al.  (54).  In  this  report,  the 
authors  assayed  expression  of  roughly  1,200  genes  in  the  same 
cell  line  after  treatment  with  significantly  higher  concentra¬ 
tions  of  lethal  toxin  and  noted  that  although  widespread  ne¬ 
crotic  cell  death  was  occurring  within  these  cultures,  essentially 
no  changes  in  gene  expression  were  detectable. 

In  contrast  to  the  trend  observed  in  unactivated  macro¬ 


phages,  we  identified  230  genes  that  changed  significantly  in 
expression  level  between  macrophages  that  were  activated  with 
LPS  and  cells  that  were  activated  and  treated  with  lethal  toxin 
(supplemental  Table  S5).  These  genes  were  sorted  by  gene 
ontology  annotation  as  described  above  in  order  to  identify 
overrepresented  functional  families  (Table  4).  In  general,  most 
of  the  host  genes  that  showed  significantly  changed  expression 
levels  following  treatment  with  lethal  toxin  were  associated 
with  signaling  and/or  control  of  transcription,  which  is  consis¬ 
tent  with  the  known  abilities  of  lethal  toxin  to  disrupt  and  alter 
host  cell  signal  transduction.  Notably,  apoptosis-related  loci 
were  not  overrepresented  in  this  set,  though  several  key  genes 
(e.g.,  ornithine  decarboxylase)  were  present,  albeit  at  lower 
levels  than  were  observed  during  infection.  We  did  observe 
some  ceil  death  in  lethal  toxin-  or  LPS-treated  macrophages 
within  the  experimental  time  course,  but  it  was  not  to  the  same 
levels  that  were  observed  in  infected  cells  (as  expected  based 
on  our  use  of  a  significant  but  sublethal  dose).  It  is  therefore 
difficult  to  assess  whether  the  absence  of  induced  apoptosis- 
related  genes  simply  reflects  lower  levels  of  cell  death  or  if 
other,  as  yet  unidentified  B .  anthracis  factors  may  play  roles  in 
inducing  apoptosis,  and  further  experiments  will  be  necessary 
to  answer  these  questions. 

Perhaps  the  most  striking  finding  observed  in  lethal  toxin- 
treated  cells  was  that  many  (142,  or  62%)  of  the  genes  that 
changed  significantly  in  activated  macrophages  treated  with 
lethal  toxin  were  also  significantly  changed  when  comparing 


TABLE  4.  Functional  families  overrepresented  in  the  transcriptional  response  to  treatment  with  anthrax  lethal  toxin  in 

activated  macrophages 


Induced 

Repressed 

Biological  process 

P 

Biological  process 

P 

Immune  response 

5.0  E-5 

Innate  immune  response  (inflammatory  response) 

0.0033 

Chromatin  assembly/modification 

1.0  E-4 

Receptor-linked  signal  transduction 

0.0043 

Regulation  of  cell  cycle 

1.3  E-4 

Regulation  of  transcription 

0.0374 

Cell  matrix  adhesion 

2.4  E-4 

Receptor-linked  signal  transduction 

5.2  E-4 

Regulation  of  transcription 

0.0149 

Intracellular  signaling  cascade 

0.0242 

Cytoskeleton  organization 

0.0253 
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FIG.  4.  Dampening  of  immune  response  genes  following  treatment  with  anthrax  lethal  toxin.  Bars  indicate  the  level  of  induction  relative  to 
expression  in  untreated  cells,  as  measured  by  DNA  microarray.  For  each  gene  indicated,  induction  levels  are  shown  for  cells  treated  with  LPS  and 
cells  treated  with  LPS  and  lethal  toxin  (LeTx). 


LPS-treated  cells  and  untreated  cells  (that  is,  they  were  also 
part  of  the  LPS-triggered  response).  Given  that  treatment  with 
LPS  appears  to  trigger  a  widespread  activation  response,  the 
significant  overlap  between  these  two  data  sets  (genes  present 
in  both  of  these  two  sets  are  noted  in  supplemental  Table  S5) 
implies  that  one  of  the  primary  effects  of  lethal  toxin  within  the 
macrophage  is  in  modulating  the  stereotypical  activation  pro¬ 
gram  induced  by  bacterial  antigens  (e.g.,  the  patterns  reported 
in  reference  6).  In  general,  many  of  the  significant  lethal  toxin- 
induced  changes  in  gene  expression  appear  to  be  in  the  down¬ 
ward  direction;  that  is,  the  induction  of  many  of  the  classical 
immune  response  genes  that  was  observed  in  LPS-treated  cells 
is  dampened  in  LPS-  and  lethal  toxin-treated  cells.  Prominent 
examples  of  this  phenomenon  are  shown  in  Fig.  4.  Previous 
studies  have  shown  that  lethal  toxin  is  able  to  suppress  a  variety 
of  signals  associated  with  activation  of  the  immune  system  (2, 
9,  10,  18,  42,  55).  Our  results  are  consistent  with  those  reports, 
and  they  demonstrate  this  effect  at  the  transcriptional  level  for 
the  first  time. 

Although  many  of  the  classical  immune  response  genes  were 
repressed  in  response  to  lethal  toxin,  a  smaller  subset  of  genes 
were  upregulated  in  lethal  toxin-treated  activated  macro¬ 
phages.  These  genes,  including  the  class  I  major  histocompat¬ 
ibility  complex  locus,  JAK-2,  calcipressin,  and  macrophage  col¬ 
ony-stimulating  factor,  are  all  associated  with  gamma 
interferon  (IFN-y)-mediated  signaling  and/or  generation  of 
T,,-l  immunity  (11,  39.  49).  IFN-y  has  been  shown  to  play  an 
important  role  in  host  responses  to  other  intracellular  bacterial 
pathogens,  including  Listeria  monocytogenes  and  Mycobacte¬ 
rium  tuberculosis  (52),  and  a  recent  study  by  Gold  et  al.  showed 
that  exogenous  IFN-7  is  able  to  improve  survival  rates  of  hu¬ 
man  macrophages  infected  with  B.  anthracis  (22).  This  report 
also  showed  that  during  infection  by  B.  anthracis ,  IFN-a/p 
signaling  is  disrupted  while  IFN-7  is  not,  which  is  consistent 
with  the  implication  from  our  data  that  although  lethal  toxin  is 
able  to  dampen  much  of  the  stereotypical  macrophage  re¬ 
sponse,  it  appears  to  leave  IFN-7  signaling  intact. 

It  seems  intuitive  that  IFN-7  and  its  associated  signaling 


pathways  would  be  important  components  in  the  macrophage 
response  to  B.  anthracis ,  and  it  is  somewhat  interesting  that 
lethal  toxin  has  not  evolved  a  means  of  IFN-y  signaling  as  it 
disrupts  other  pathways.  Although  IFN-y  has  a  protective  ef¬ 
fect  on  infected  macrophages  in  cell  culture  and  thus  appears 
to  work  to  the  detriment  of  the  infecting  bacteria,  it  remains 
possible  that  there  are  more  broad-reaching  effects  of  IFN-y 
signaling  in  vivo  that  may  work  to  the  advantage  of  B .  anthracis 
at  the  organismal  level  and  that  these  effects  ensure  that  there 
is  no  selective  pressure  toward  evolution  of  a  means  for  IFN-y 
disruption. 

Conclusions.  This  study  represents  the  first  characterization 
of  the  transcriptional  changes  within  macrophages  during  in¬ 
fection  by  Bacillus  anthracis  and  the  first  genomewide  charac¬ 
terization  of  transcriptional  changes  in  macrophages  treated 
with  anthrax  lethal  toxin.  Although  many  of  the  genes  identi¬ 
fied  as  differentially  regulated  have  been  previously  shown  to 
play  a  role  in  the  general  innate  immune  response,  we  have 
also  identified  a  limited  set  of  macrophage  genes  that  are 
induced  or  repressed  in  response  to  B.  anthracis  but  have  not 
yet  been  identified  as  such  in  other  bacteria-host  systems. 
These  genes  provide  a  transcriptional  basis  for  many  of  the 
phenomena  observed  during  anthrax,  including  the  induction 
of  apoptosis  within  host  phagocytes,  and  it  will  be  of  consid¬ 
erable  interest  to  see  if  future  studies  characterizing  the  mac¬ 
rophage  response  to  other  pathogens  show  the  regulation  of 
any  of  these  genes  to  be  truly  specific  to  the  B.  anthracis- 
macrophage  interaction. 

We  anticipate  that  this  analysis  of  gene  expression  patterns 
during  infection  and  toxin  treatment  will  be  valuable  in  ongo¬ 
ing  efforts  to  develop  new  detection  and  therapeutic  options 
for  anthrax.  It  must  be  noted  that  these  studies  were  all  done 
with  an  attenuated  strain  of  B.  anthracis ,  and  thus,  although  the 
data  presented  here  are  comprehensive  relative  to  that  strain, 
fully  virulent  B.  anthracis  (carrying  pX02)  may  cause  other 
changes  in  host  gene  expression  that  were  not  identified  here. 
It  will  therefore  be  of  considerable  interest  to  extend  these 
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studies  to  wild-type  B.  anthracis ,  and  we  are  pursuing  these 
studies. 

We  expect  that  these  data  will  be  of  use  to  those  in  the  B. 
anthracis  field  as  well  as  to  the  general  pathogenesis  commu¬ 
nity,  and  we  have  made  all  microarray  data  collected  in  this 
study  publicly  available  via  the  TIGR  website  (http://pga.tigr 
.org)  and  in  the  ArrayExpress  database  (http://www.ebi.ac.uk 
/arrayexpress/,  accession  number  E-TIGR-84).  We  are  con¬ 
tinuing  to  examine  these  data  for  important  insights,  and  we 
are  also  extending  the  techniques  described  in  this  work  to  the 
study  of  primary  cells  in  culture  and  in  vivo. 
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